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1. Introduction 

Many precision experiments are increasingly becoming sensitive to nonperturbative QCD cor- 
rections. For example, the measurement of the magnetic moment of the muon currently shows a 
discrepancy with the standard model of over three standard deviations. Its theoretical uncertainty is 
dominated by hadronic effects, making the muon g - 2 a prominent example of the importance of a 
fully nonperturbative determination of hadronic corrections. This is just one example. The signif- 
icance of QCD corrections to otherwise precision observables will increase with the experimental 
programs envisioned for the future. In many cases, the discovery of physics beyond the standard 
model may depend on accurate control of these hadronic effects. 

In these proceedings, we will discuss several opportunities for lattice QCD calculations of 
hadronic corrections to important measurements that may in fact be more feasible than previously 
thought. The QCD corrections to the muon g-2 will serve as a concrete example. This will allow us 
to identify an issue that makes the calculation of these quantities more difficult and then describe a 
modified method that was introduced to alleviate this problem [1]. Additionally, there is a growing 
lattice effort on precisely this quantity and we will review the latest results. 

After having laid the groundwork with the muon g-2, we will continue with several addi- 
tional examples that illustrate the potential impact of the modified technique. We will examine 
the leading-order hadronic corrections for the electron and tau magnetic moments and the leading 
QCD contributions to the running of the QED coupling. Additionally, we will note applications to 
the Adler function, the determination of the strong coupling constant and the QCD corrections to 
the energy levels of muonic-hydrogen. We will then close with a calculation of the next-to-leading- 
order vacuum-polarization corrections to the muon g-2. 

2. Leading-order QCD correction to the muon magnetic moment 

The BNL measurement of the anomalous magnetic moment of the muon = (g^ — 2)/2 [2] 
and the standard model estimate thereof [3] differ by more than three standard deviations. This 
discrepancy may indicate physics beyond the standard model, but making such a statement defini- 
tively requires a thorough understanding of all sources of uncertainty and ideally a significantly 
larger discrepancy. The experimental community is pursuing two future muon g-2 experiments at 
Fermilab [4] and J-PARC [5], aiming to improve the experimental precision on from 6.3 • 10 
to (1—2) -10~ 10 . Since 1.2- 10~ 3 , the new experiments will reduce the relative precision from 
0.5 • 10~ 6 to 0.9 • 10~ 7 . At this precision, the comparison between theory and experiment would be 
dominated by the standard model uncertainties alone, hence improvement from the theory side is 
highly desirable. 

The value of receives contributions from all parts of the standard model, each contributing 
to the theoretical uncertainty as shown in table 1 . Quite clearly, the standard model uncertainty is 
overwhelmingly dominated by hadronic physics. The total QCD contribution <3^ CD can be orga- 
nized as an expansion in the QED coupling a as follows 

^ CD = a 2 Af + a 3 4 3) + ^(« 4 ) = af + af + ^(a 4 ) , (2. 1) 
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Contribution 


Error [1(T 10 ] 


QCD-LO 


5.3 


QCD-NLO 


3.9 


QED/EW 


0.2 


Total 


6.6 



Table 1: Standard model uncertainties in [3]. 
The contributions QCD-LO and QCD-NLO refer 

(2) (3) 

to fljj and ah_ in equation 2.1. All remaining 
contributions are collected together and labeled 
QED/EW. 



Y 




Figure 1: Leading-order QCD contribution ah . 

(2) 

The QCD contribution to a„ , denoted by the 
shaded region, can be related to R{s), given in 
equation 2.2, determined experimentally or calcu- 
lated using lattice QCD. 



where = (XML . The expansion in a is perturbative but the must be calculated nonpertur- 

(2) 

batively. The largest source of uncertainty is due to ah , which we discuss next. 

(2) 

2.1 Experimental determination of a^ 

(2) 

The leading-order correction ah , shown in figure 1 , can be determined experimentally. It can 
be written as an integral of R(s) and a known function K(s/mj l ) as [3] 

(2) ? r° ds , , , , , , , , a(y* — > hadrons) 
<A 2) = a 2 / —K(sml)R(s) with R(s = — . (2.2) 

In practice, measurements of a(y* — > hadrons) from many different experiments are combined to 
form the integral above. The most recent compilation of measurements results in = 6.923 (42) • 
10~ 8 [6], which is an improvement on the error given in table 1. This approach has been and will 
continue to be for some time the most accurate means of providing the QCD input needed to form 
the standard model prediction for a^. This value is significantly more precise than current lattice 
results, but we should bear in mind that this result requires a substantial experimental effort to 
determine a quantity that should in principle be predicted from the theory itself. Reaching and 
even exceeding the precision of the e> 
term efforts of the lattice community. 



(2) 

even exceeding the precision of the experimental determination of ah should be part of the long 



(2) 

2.1.1 Estimates of the flavor dependence of a^ 

As we will see shortly, lattice calculations must include the four lightest quarks to reach the 

(2) 

precision on needed for the future muon g-2 experiments. However, the number of quark 
flavors used in current lattice computations varies from «y = 2 to 3 and just recently 4 flavors. 

(2) 

Thus it is useful to have some guidance on the rif dependence of ah . Unfortunately, there is no 

(2) 

unique way to do this for ah , which is a problem for many other observables as well. The effects of 
decoupling a quark flavor depend on the renormalization conditions used for the remaining degrees 
of freedom. This ambiguity exists in perturbation theory and would also apply to any analysis of the 
experimental results. The advantage of lattice calculations is that they can provide a well-defined 
way to prescribe such a definition. 
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E [GeV] 

Figure 2: Measured R(E). The R(s) with y/s = E 
from F. Jegerlehner [7] is shown. 



tlf 


a [ ^n f [Exp] 


a% f [Lat] 


5 


6.93 (06) 




4 


6.93 (06) 


underway [8] 


3 


6.81(05) 


6.41(46) [9], 6.18 (64) [10] 


2 


5.67 (05) 


5.72(16) [1], 5.46(66) [10] 



Table 2: Estimated n / dependence of (given 
in units of 10~ 8 ) from experimental measurements 
using equation 2.3 ("[Exp]") and current lattice cal- 
culations ("[Lat]"). 



Keeping in mind these limitations, we proceed with a simple means of estimating the flavor 
contributions from the experimental measurements. We start with the experimentally determined 
R(s) [7] and rescale by the electric charges Qf of the relevant quark flavors / as follows 

R nf (s) = R{s) -t—f for W n < * < with a% f = a 2 f —K(s/m 2 )R nf (s) . (2.3) 

The YJ} accounts for the n degrees of freedom present in the experimental measurements between 
the quark thresholds Am\ and 4ra^ +1 and the Y!7 restores the desired rif flavors. The leading-order 
perturbative contribution to R(s) scales this way, so this prescription is valid up to perturbative 
corrections and away from the resonance regions. Near resonances, this amounts to a quark-hadron 
duality argument. For example, the strange quark content of the prominent meson would not 
follow such a scaling but the integrand K(s/mj l ) is relatively smooth and effectively averages R(s) 
across a window in s over which this scaling is expected to be more effective. In practice, the part 
of this prescription that has the most impact on the resulting estimates of '„ is the use of quark 
masses rather than a resonance mass to define the thresholds. And in fact, only the strange quark 
threshold is particularly sensitive to this choice. 

(2) 

Using this simple prescription, the rif dependence of can be estimated. To do this, we 
have used the R(s) compiled by F. Jegerlehner [7]. This is shown in figure 2. The results of 
this procedure and a comparison with the current lattice calculations are given in table 2. This 
suggests that the charm quark correction is roughly 1.2 • 10~ 9 . The current experimental precision 
is 6.3 • 10~ 10 , which makes it clear that charm quark contributions are already necessary to reach 
the precision of the BNL measurement, let alone the precision for the future g-2 experiments. 

Alternatively, one could attempt to assign flavor weights to each of the final states in the total 

(2) 

cross section a(y* — > hadrons) and form a^' nf . This approach is an alternative prescription to the 
one given above and, in fact, differs from it [11, 12], illustrating the ambiguity in extracting the 
rif = 2 piece. It turns out to lead to a larger rif = 2 contribution than given in recent nj = 2 lattice 
calculations, but we want to emphasize that any such ambiguities are systematically eliminated as 
the lattice calculations account for all the relevant quark flavors, which appears to be rif = 4. 

(2) 

2.2 Lattice calculation of a^ 

(2) 

The standard method to calculate ah using lattice QCD was given by Blum in [13]. It requires 
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Figure 3: Vector-meson coupling g V - 
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Figure 4: Modified method a 



(2) 



calculating the vacuum-polarization function H(Q 2 ) and evaluating the integral 

dQ 2 



a 



o Q 



2 co{Q 2 /ml)U R {Q 2 ) 



(2.4) 



The weight function CO accounts for the perturbative portion of the diagram in figure 1 and is known. 
ns(<2 2 ) is the once-subtracted vacuum-polarization function, Hr(Q 2 ) = n(<2 2 ) — 11(0), where n 
is given by 

(Q,Qv-Q 2 d, v )U(Q 2 )^ f d'Xe'^-^iJ^MY)) , 



and is directly calculable in Euclidean space. This formulation can be related to the approach used 
to experimentally determine ai by noting that R(s) is proportional to ImII(— s + ie), which is 
non-zero only on the branch cut starting at s = 4m 2 . A standard dispersion analysis then allows 
one to relate equations 2.2 and 2.4. 

(2) 

2.2.1 Role of external leptonic scale in 

(2) 

Naively, the lattice calculation of should be relatively easy. The only non-trivial part of 
the computation is the calculation of the Euclidean two-point correlation function (J I1 (X)J V (Y)}, 

(2) 

which can be accurately determined. Furthermore, the quantity is dimensionless, so it seems 
reasonable to expect that it may exhibit a relatively mild dependence on the scales in the problem, 
particularly the quark masses and lattice spacing. (Figure 3 shows an example of a typical dimen- 

(2) 

sionless quantity.) However, is dimensionless only at the expense of introducing an external 
scale, the muon mass m^, and this has several consequences. 

First, win introduces a dependence on the lattice spacing a in an otherwise dimensionless ob- 

(2) 

servable. We can see this by writing the integration variable for in lattice units 



a 



dO 2 

^ w (G 2 /(^) 2 )n la t(<2 2 ) 



where Q = aQ is the momentum variable in lattice units and rii at (<2 2 ) = Hr(Q 2 /a 2 ) is directly 
calculated in lattice units. Thus the lattice spacing a is needed in physical units to form am^. This 
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suggests that a„ may behave more like a dimensionful quantity. As a second consequence, the 
introduction of m^ also allows for a stronger quark mass dependence than might otherwise be ex- 
pected. The dominant contribution of the lightest vector-meson with mass my and electromagnetic 
coupling gv is proportional to gymtfrriy. This is a model-dependent statement, but it is suggestive 
that may in fact behave more like a mass dimension -2 observable. (A plot of my corresponding 
to the gy in figure 3 is given in [1].) 

These two observations can be made precise by introducing an effective dimension 

a dX 

This quantity is defined so that it is sensitive to only QCD scales rather than the overall dimension. 
To accomplish this, we write the observable X as a function of both the lattice spacing a and the 
coupling g separately X =X(a,g). Of course, a = a(g) is eventually chosen to be some function of 
the coupling, but treating a and g separately allows us to isolate the impact of the scale setting on 
the quantity X. Furthermore, fi? e ff is defined so that it reproduces the usual definition of dimension 
for a simple QCD observable but it differs for composite observables. 

Several examples may help illustrate c/ e ff . First consider some observable M that is a standard 
QCD quantity of mass dimension n. This quantity would satisfy M(a,g) = a~ n M(g), where M(g) 
is what is calculated on the lattice and the factor of a~" is eventually put in by hand. Then for such 
an M we have 

Thus the dimension of quantities that we normally calculate is unaltered. However, for a composite 
object the answer can differ. For example, 

d e f[ [gv m l/ m v] = d eS [gva 2 m* /m v ] = -2 

where gyffti/my is the leading piece of the vector-meson contribution to mentioned earlier. 
Thus d^ft captures the dimensionality of the QCD scales in this rather simple composite expression. 

(2) 

We can now apply the definition of d^s to ah and we find 

d^[af] = -2^[ d ^co(Q 2 /ml)Q 2d ^ ^*{&/n$U*{&)} . 

(2) 

This quantity has a continuum limit. It is rather easy to show that <5f e ff [a\i ] < 0, making it clear that 

(2) 

Ofi acts like an observable with a negative mass dimension. Additionally, for — > 0, we have 
<sf e ff — > —2, and for m^ — > °°, we can show that d e f[ — > 0. For an intermediate mass, this quantity 
must be calculated nonperturbatively. For the muon, we find 

4ff[«t 2) ] = -l-887(5) 

(2) 

which clearly indicates that behaves much more like a mass dimension -2 quantity than a 
dimensionless observable. 
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(21 

2.3 Modified lattice method for ; 

(2) 

Having diagnosed the difficulty in , first with dimensional analysis and then a model ar- 
gument, and having provided a clean definition of the problem using d e &, we can now attempt to 

(2) 

remedy it. In the end, we will define a modified quantity a\y that has the same physical limit as 

yet satisfies d e ff [aL ] = 0. Since the physical values of both observables are the same, we can 

safely use either quantity to perform the computation. Furthermore, our observation that J e ff = 

(21 

will provide a theoretical explanation for why ay should lend itself to an easier calculation. 

(21 (2) 

Starting with the observation that c/ e ff [a\i ] ^ 0, we sought a minimal way to modify a^ to 
eliminate this unexpected dependence on the lattice spacing. This is caused by the fact that is 
an external scale and not capable of absorbing the dependence on a. The solution we came upon 
was to insert the factor H^ hys /H 2 as follows 



in 



A' 



where H is some hadronic scale and // p h ys is its physical limit value in physical units. Additionally, 

(2) 

the value of H is understood as being calculated at the same mp$ as . Each choice of H gives 

12) 

rise to a distinct new observable ay. We will shortly pick a favored H, so the dependence on H 

12) 

in defining a.y is suppressed. By construction, this quantity has the same physical limit as the 
standard definition and eliminates the unwanted dependence on the lattice spacing, 

(2) (2) (2) 

lim a n —ah' and d e nWn-) = 0. 

This modification completely eliminates the unwanted a dependence, but it does not automatically 
weaken the quark mass dependence. We know that the vector-mesons make a dominant contribu- 
tion to aj][ of the form g\m 1 ^jm\ and that gy, shown in figure 3, is only weakly nips dependent. 
This suggests choosing H = my. Other choices have been examined in [1] but we will use H = my 
exclusively in these proceedings when discussing the modified approach. 

(2) (2) 

The results for a-^ are given in figure 4. As a,y is now a proper dimensionless observable 
composed of only QCD scales, it behaves like any other dimensionless quantity, as we can check 
by comparing with the gy in figure 3. Additionally, the choice H = my absorbs much of the mps 
dependence that has troubled previous calculations using the standard method. All indications are 
that the systematic errors are relatively mild. Additionally, we have checked that disconnected 
diagrams do not rise above the statistical errors shown in figure 4. The extrapolated value is given 
in table 2 and is consistent with the estimated «/ = 2 piece of the experimental determination. More 

(21 

importantly, it is encouraging that the resulting error on the physical limit value of a^ is already 
within a factor of 3—4 of the experimental determination using equation 2.3. 

(2) 

2.4 Comparison of current lattice calculations of a^ 

(2) 

Several groups have performed calculations of a^ with both nj = 2 [1, 10] and rif = 3 [14, 9], 
and a first tif = 4 calculation is underway [8]. To compare results, we focus on the standard method 
for this section. First, we can examine the lattice calculations for n/ = 2, which are shown in fig- 
ure 5. There we find agreement between both computations, including their extrapolated values, 
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O Imp. Clover, N f =2 



0.1 
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0.3 



0.4 



(2) 

Figure 5: Two-flavor lattice calculations of . 
The calculations are from [1] ("Twisted") and 
[10] ("Imp. Clover"). 
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Figure 6: Three-flavor lattice calculations of 
4 2) - The results are from [14] ("Asqtad"), [9] 
("DWF") and [10] ("Imp. Clover"). 



which are given in table 2. These results were calculated using different actions, at least two lattice 
spacings, multiple physical volumes, a broad range of pion masses, and different treatments of the 
low Q 2 extrapolations of H(Q 2 ). The level of agreement is rather compelling for these two calcu- 
lations. Additionally, it appears that the calculation of [10] shows a rapid rise as mp$ is lowered, 
which is consistent with the model-dependent expectations from the vector-meson contribution 
discussed previously. 

We now turn our attention to the rif = 3 calculations shown in figure 6. The situation here 
appears to be less compelling than for rif = 2, but we must take some care before reaching such 
a conclusion. A detailed comparison requires a bit more space than allowed in these proceedings, 
so we will limit ourselves to commenting on one conceptual issue that is important to understand 
before making any definitive statements. When comparing the rif = 3 results, we must consider 
how the strange quark mass is determined. Unless the chosen renormalization conditions match, 
there is no reason why these curves would generally agree. The only expectation is that the values 
extrapolated to nips = m n must agree when all other uncertainties have also been controlled for. 
This seems to be the case in figure 6. Additionally, the extrapolated values given in [9] and [10] 
(also given in table 2) do agree. (The work of [14] did not cite a final result but it overlaps with 
the calculation of [9] for all three values of nips used in [14].) Further understanding is needed for 
rif = 3, but the results of the lattice calculations may be more encouraging than is reflected in a 
simple head-to-head comparison as done in figure 6. 

3. Leading-order QCD correction to the electron magnetic moment 

We now turn our discussion to a sequence of examples illustrating the application of the mod- 
ified method of [1]. The first two examples are simple extensions to the other two leptons but 
nonetheless provide nontrivial tests of the new method. Besides, these are first lattice calculations 
of both quantities. In [1], the leading-order QCD contributions to the electron and tau magnetic mo- 

(2) (2) 

ments were also calculated, a e and a\ . The electron magnetic moment is measured to a precision 
of 0.28 parts-per-trillion [15]. The lightness of the electron makes a e significantly less dependent 
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Figure 7: Modified method a 



(2) 



Figure 8: Modified method AE. 



on QCD corrections but the enhanced precision of the experimental measurements overcomes the 
reduced sensitivity. The current measurement of a e is so precise that it is in fact used to determine 
a. The standard model is then tested by comparing this value to other determinations of a. These 
comparisons have reached the precision that QCD effects of a e can not be ignored. However, the 

(2) 

precision on a e is not yet high enough to probe the error on a e , so there is no pressing experimental 

(2) 

need for higher precision a e determinations. 

(2) 

For our puiposes, a e provides a test of the modified method that is sensitive to only the 
extreme lowest Q scales. To a high accuracy, a e can be approximated as 



3 e dQ 2 



+ tff(m e /A) A 

=0 



(2) 

This approximation is not used in [1], but we mention it in order to illustrate that a e probes 
essentially just the derivative of H(Q 2 ) at Q 2 = 0. It also illustrates in a simple way again the 

(2) (2) 

idea that a e may behave differently from a typical dimensionless observable. In the case of a e , 
the approximation above strongly suggests that it will behave very much like a mass dimension -2 
observable. This can be made precise by noting that 

<4ff[4 2) ] = -1.999984(1). 

(2) (2) 

The results for a e from [1] are shown in figure 7 using the modified method. The behavior of a e 

(2) 

is very similar to and agrees with the estimated «/ = 2 piece of the experimental measurement. 

Additionally, the leading QCD vacuum-polarization correction to the 2P—2S Lamb shift of 
muonic-hydrogen is also proportional to the slope of n(<2 2 ) at Q 2 = and hence closely related to 

(2) 

a e . This correction is given in [16], with m r the reduced mass of the jl—p system, as 



AE = 27ta 5 ml dIl « 
' dQ 2 



Q 2 =0 



3 dn R (Q 2 /H 2 hys .H 2 

and AE = Inarm* ^ 

dQ 2 



Q 2 =0 



The second form is the modified approach for AE and results from a consistent treatment of external 
scales as discussed thoroughly in section 5. The results for AE are shown in figure 8. 
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Figure 9: Modified method a 
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Figure 10: Modified method Aa (1) (M^). 



4. Leading-order QCD correction to the tau magnetic moment 



The magnetic moment of the tau is substantially more sensitive to potential new physics than 
that of the muon, but it has not been experimentally measured. There are experimental bounds on 

(2) (2) 

a % [17], but more interesting to us is that a\ can be determined by the same analysis used for 

(2) 

and a e . Additionally, due to the heaviness of the tau, it is sensitive to a different range of QCD 
scales. In fact, 



( 2 )l 



-0.936(13) 



(21 (2) 

indicates that it behaves quite a bit different than and a e and more like a mass dimension 

(2) 

-1 observable. In fact, applying the standard method to a\ leads to a reasonable calculation. 
However, the arguments for the modified method still apply here and it is yet another test of the 

(2) 

approach. The results for a\ are shown in figure 9. Again, the resulting observable has a mild nips 
dependence. Furthermore, it agrees with the estimated «/ = 2 piece of the dispersive result, deter- 
mined from equation 2.3, providing more confidence in both the modified method and the simple 
prescription for analyzing the tif dependence of R(s). As figure 9 shows, the lattice determination 

12) 

of a\ is already more accurate than its experimental determination, suggesting that future mea- 

(2) 

surements of a x can in fact rely on a fully nonperturbative determination of a\ without sacrificing 
any precision once rif = 4 lattice calculations have been completed. 



5. Leading-order QCD contribution to the running of the QED coupling 

In order to further demonstrate and understand the modified method of [1], we have applied 
the same idea to the determination of the leading-order QCD corrections to the running of a. This 
is the first lattice calculation of this quantity and all results discussed in these proceedings are 
preliminary. The running of a is normally treated by introducing an effective coupling given by 
summing all one-particle irreducible bubble insertions in the photon propagator. This results in 
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The QCD contribution AaQ CD is again expanded in a. The leading-order correction is 

Aa^(Q 2 )=47uaU R (Q 2 ). 

The value of a(Q 2 = 0) is just the usual coupling a, which is known to a relative precision of 
1 • 10~ 9 . However, after evolving a to a high scale, say the Z-boson pole at Q 2 = Aff, the relative 
precision on ce(M|) drops to 1 • 10~ 4 [6], making Of(M|) one of the more poorly known fundamen- 
tal parameters in high energy predictions. Similar to an, the dominant uncertainty in this evolution 
is due to hadronic corrections, which is then passed into every high energy process through the use 
of the running coupling a(M|). This has larger impact than one might have naively expected. As 
one example, a recent global analysis by the Gfitter collaboration determined the Higgs mass to 
be ntH = 44I43 GeV, if the experimental determination of Ad£ (M§) was not included in the fit and 
found m H = 96^ GeV if it was [18]. 

(2) 

The treatment of the external scales for a) introduced earlier uniquely fixes the treatment of 
the now external scale Q 2 in Aa(Q 2 ). To see this, we can simply rewrite equation 2.5 as 

af = a 2 jf d -^co (Q 2 /m 2 ) U R (Q 2 /H 2 hys • H 2 ) . 

This suggests rather clearly that we should consider the following modified definition of Act 1 ) 

A«( 1 )(2 2 ) = 47ian R (Q 2 /H 2 hy ,-H 2 ) . 

(2) 

Just as for cij , this new quantity explicitly has the correct physical limit but also satisfies d e s = 0. 

We can examine the consequences of this definition by first focusing on Q 2 = Mq with Mq = 
2.5 GeV, which is a common matching scale for phenomenological work. The lattice calculation 

(2) 

is shown in figure 10. We see that, just as for each of the a\ calculations, the modified def- 
inition results in a rather mild looking extrapolation to the physical point, giving Aa^iM^) = 
5.72(12) • 10~ 3 . We can also apply the same treatment of the tif dependence of the experimen- 
tally determined R(s) to AaW resulting in Aa (1 )(M ( 2 ) = 5.60(06) • 10~ 3 . The preliminary lattice 
computation results in an uncertainty that is now only twice the error of the experimentally deter- 
mined quantity, suggesting that lattice calculations could be a competitive, if not superior, way to 
determine Aa(Q 2 ) at least for low scales. 

Now, we can repeat the analysis for all Q 2 and determine the QCD induced running of a(G 2 ). 
This is shown in figure 1 1. By comparing results with different lattice spacings, we notice signifi- 
cant lattice artifacts only for Q 2 > 7 GeV 2 . This appears to be just a mild obstacle to an accurate 
determination of Aa(<2 2 ) in the relevant low Q 2 regime. The scale of Q 2 = M 2 was chosen be- 
cause perturbation theory becomes reliable for yet larger Q 2 . To run a to higher scales, we have 
determined a s by matching Yl{Q 2 ) to the perturbative expectations for the Q 2 regions that can be 
reached in lattice calculations and then determined Aa^'(2 2 ) at larger scales through 

Aa (1) (M|) = Aa (1) (M ( 2 ) + (Aa (1) (Ml)-Aa( 1 )(M^)), 

where the perturbative expression for Aw 1 ' (M|) — Ao^ 1 ) (Mq ) is available at 5 loops [19]. 
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We note that to determine a s , it appears to be better to calculate the Adler function D(Q 2 ) = 
dYl{Q 2 ) /d\n(Q 2 ). To consistently apply our treatment of the external scale Q, we define a modified 
Adler function 

D(Q 2 )=D(Q 2 /rf hys -H 2 ). 
The results of a preliminary calculation with the modified technique are shown in figure 12. 

6. Next-to-leading-order QCD contribution to the muon magnetic moment 

The precision of the current BNL measurement of already requires that the next-to-leading- 
order QCD correction be accounted for. Most of the diagrams involve insertions of the vacuum- 
polarization correction into lower order QED diagrams, but a new QCD contribution, called light- 
by-light, also occurs at this order. 

6.1 Vacuum-polarization corrections 

The vacuum-polarization correction can be inserted once or twice into any photon line of a 
two-loop or one-loop QED diagram, respectively. This results in 16 diagrams that involve one oc- 
currence of n(<2 2 ) and one diagram with two insertions of Yl(Q 2 ). Example diagrams are shown 
in figure 13. Expressions in various forms are available for these corrections when expressed as 
integrals over R. The analytic continuation to Euclidean space is complicated for some of these 
contributions, but it appears that all vacuum-polarization contributions to can be calculated in 
Euclidean space. A preliminary calculation using the modified approach 1 , shown in figure 14, gives 
fl (3,vp) _ _7 99 . 10~ 10 . The uncertainty includes only statistical errors but the systematic un- 
certainties appear to be small. For comparison, the tif = 2 piece of the experimental measurement 
is a^ 3 ' vp ' = —7.94(16) • 10~ 10 . Further study is underway, but the initial results for the vacuum- 
polarization contribution to seem to agree with the expectations from the experimental mea- 
surements and have a nearly comparable uncertainty. Most importantly, the precision of the lattice 
result is better than the accuracy expected for the future muon g-2 measurements, so it seems that 

'At the conference, only a partial accounting of the diagrams was given and resulted in a different value. The result 

(3) 

reported here accounts for all vacuum-polarization contributions to . 
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Figure 13: Example diagrams for 
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lattice QCD should be quite capable of determining the higher-order vacuum-polarization contri- 
butions at the required precision. 



6.2 Light-by-light corrections 



(3) 

In contrast to the vacuum-polarization corrections, the light-by-light contributions to au rep- 
resent a real challenge. There are ongoing lattice studies by several groups [20, 21, 22], each using 
different methods but still all exploratory. The uncertainty on the light-by-light correction a^' lbl ^ is 
nearly as large as that of and there are open questions regarding the methods currently used for 
its determination, so a nonperturbative calculation of a^' M ^ is highly desirable. However, a lattice 
calculation of one piece of the higher-order contribution is less satisfying than a complete lattice 
calculation of both the leading-order and next-to-leading-order corrections. It now appears that the 
vacuum-polarization pieces should be calculable, so the light-by-light contribution is the only re- 
maining piece needed for a completely nonperturbative determination of <3^ CD accurate to &{a i ). 
We can only hope that this will encourage an even greater effort within the lattice community to 
tackle the desperately needed light-by-light contribution. 



7. Conclusions 



We have discussed several examples of important measurements that receive sizable hadronic 
corrections. The muon g-2, which hints at beyond-the-standard-model physics, is the most press- 
ing observable, but Aa(Q 2 ) may in fact have a much broader impact on precision standard model 
predictions. Using a modified lattice approach, the leading QCD corrections to both of these ob- 
servables appear to be reliably calculable. To further explore the new method, we have also exam- 
ined the leading corrections to the electron and tau leptons and the Lamb shift in muonic-hydrogen. 
We have calculated the Adler function, which can be matched to perturbation theory to determine 
the strong coupling. Lastly, we have worked out methods to examine all the vacuum-polarization 
corrections at the next-to-leading order. In several cases, the currently reached precisions on these 
quantities are approaching that of the corresponding experimental determinations. This indicates 
that fully nonperturbative determinations of QCD corrections to electroweak observables may be 
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feasible at the precisions needed by future experimental measurements that aim to discover or con- 
strain physics beyond the standard model. 
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